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Abstract 



We search for direct CP violation in flavor specific B Kir decays by mea- 
suring the rate asymmetry between charge conjugate modes. The search is 
performed on a data sample of 11.1 million BB events recorded on the T(4S') 
resonance by the Belle experiment at KEKB. We measure 90% confidence 
intervals in the partial rate asymmetry ^cp of —0.25 < A(jp{K^tt^) < 0.37, 
-0.40 < AcpiK^TT^) < 0.36, and -0.53 < ^cp(i^°vr=F) < 0.82. By 
combining the K^-ir^ and K^tt^ final states, we conclude that —0.22 < 
-4cp(i^'F(7r=^ + 7r°)) < 0.25 at the 90% confidence level. 
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Typeset using REVT^ 



The most straightforward indication for CP violation in the B meson system would be 
a time-independent rate asymmetry between CP conjugate decays into flavor specific or self 
tagging final states. Direct CP violation (DCPV) [|^ of this type will occur in a decay 
containing at least two amplitudes that have different CP conserving and CP violating 
phases. The charge asymmetry will be most evident in decays where the two amplitudes are 
of comparable strength. The partial rate asymmetry can be written as 

A, 
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+ N{B - 




- + 2IA1IIA2I cos5cos0 



CP 

where S and are the CP conserving and CP violating relative phases between amplitudes 
Ai and A2] B represents either a B^ or B'^ meson, / represents a self tagging final state, 
and B and / are the conjugate states. 

In the Standard Model, DCPV occurs in charmless hadronic B decay modes that involve 
both penguin (P) amplitudes and weak h ^ u tree (T) amplitudes containing the CP 
violating weak phase 03 = aig{y*^) (in a standard convention 0). Two examples are 
B^ —>■ K^TT" and B^ K^tc^. Uncertainties in the relative size of the P and T amplitudes 
as well as their relative strong phase make it difficult to predict the magnitude of Acp- 
Attempts have been made to estimate or calculate these quantities however, the 

different methods give conflicting results. Although absolute predictions vary, they tend to 
agree that ^cp has approximately the same magnitude and sign in the K'^n^ and K'^n^ 
final states This suggests that we can combine these two modes to increase the 

statistical precision. 

At the quark level, the B^ — * K^tx^ process does not include a contribution from the 
tree diagram; thus one expects a negligible asymmetry in this mode. DCPV can be induced 
in this mode through rescattering effects 0; however, these are expected to be small 

In this paper, we describe a search for DCPV by measuring the difference between the 
yields of B and B decays into the self tagging final states i^'^vr^, fT^vr^, and K^n^ . Searches 
for partial rate asymmetries in the B —>■ Kir decays have been reported in DCPV has 
been observed in the kaon system P,p!0|]. 

The analysis is based on data taken by the Belle detector [jll| at the KEKB [|l2l e~^e~ 
storage ring. The data set consists of 10.4 fb~^ on the T{AS) resonance, which corresponds 
to 11.1 X 10^ BB events. A 0.6 fb~^ data sample is taken 60 MeV below the T{AS) resonance 
to perform systematic studies of the continuum e~^e~ qq background. KEKB is a two- 
ring, energy-asjTiimetric storage ring with electrons stored at 8 GeV and positrons stored 
at 3.5 GeV giving an T(4S') center-of-mass system (cms) boosted by 7/? = 0.425 in the 
electron beam direction. The Belle detector is described in detail in |jlT|. Charged vr and K 
mesons are identified based on their energy loss [dE/dx) in the central drift chamber and 
their Cerenkov light yield in an array of silica Aerogel counters (ACC). For each K and vr 
hypothesis, the dE/dx and ACC probability density functions are combined to form a single 
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likelihood, Lk or £,r- ^ and vr mesons are distinguished by cutting on the likelihood ratio 
I^kI{,I^k + -^^vr)- The particle identification (PID) efficiency is determined using nearly pure 
samples of K and vr mesons tagged using the continuum D*^ decay chain D*^ D^tc^, 
— s> K^TT^ For tracks in the B ^ hh {h = K ot it) signal region of 2.4 GeV/c 

< Perns < 2.8 GeV/c (1.5 GeV/c < p\ah < 4.5 GeV/c) the K efficiency and fake rate are 
ex = 0.85 and J'k = 0.10 (true K fakes tt); the vr efficiency and fake rate are 6-,^ = 0.92 and 
/,r = 0.04 (true vr fakes K). The relative systematic error is 2.5% in the efficiencies and 10% 
in the fake rates. 

In the K'^ir^ modes, the asymmetry is diluted by double mis- identification of K^ti" as 
Ti^K~ . The dilution is given by 

= (1 - 2w)All\ 

where w = f / (e + f) is the wrong tag fraction, / = fKfn = 0.0040 ± 0.0006 is the double 
mis-identification rate and e = ex^n = 0.78 ± 0.03 is the PID efficiency. This gives a wrong 
tag fraction of w = 0.0051 ± 0.0007 and a dilution of 1 - 2u; = 0.990 ± 0.002. The value of 
AcpiK'^n^) is corrected by this factor. 

The identification procedure of i? ^ Ktt decays has been documented in Jill. Here we 



give a brief overview of the analysis. All charged tracks that originate from the interaction 
point are considered as kaon or pion candidates. Kg mesons are reconstructed using pairs 
of charged tracks that have an invariant mass with |m(7r"''7r~) — m{K^) \ < 30 MeV/c^. The 
candidate must have a displaced vertex and flight direction consistent with a high momentum 
Kg originating from the interaction point. Candidate vr'' mesons are formed by requiring 
pairs of 7s to have |m(77) — m(7r°)| < 16 MeV/c^. Continuum background is reduced by 
cutting on a likelihood ratio formed from event topology variables and initial state angular 
momentum variables such as the B flight direction or the decay axis direction. This cut 
removes more than 95% of the continuum background while retaining 40 to 50% of the 



signal depending on the flnal state. Further details can be found in [ITl 



To reconstruct the B meson, we form two quantities: the energy difference, AE = Eb — 



Eheam, and the beam constrained mass, m^c = \J E'^ea.m ~ V^bi where -Ebeam = = 5.29 

GeV, and Eb and pb are the reconstructed energy and momentum of the B candidate in the 
cms. For the K^tt^ and K^gir"^ ffnal states, the signal yield is extracted by fitting the AE 
distribution in the region m^c > 5.27 GeV/c^ and \ AE\ < 250 MeV. As the AE distributions 
for the K^Tx^ final states have large asymmetric tails, we apply a two-dimensional fit to the 
AE and rrihc distributions in the region m^c > 5.2 GeV/c^ and —450 MeV < AE < 150 
MeV. The Gaussian signal shapes are calibrated using the B^ D^-k^ data sample and high 
momentum K^tt^, K'tc'^tt^, and —>■ K^tc'^ decays where the daughter particles 

have the same momenta as signal particles [jl3|. The continuum background shapes are 



determined from sideband data. We use a linear function for the AE background; we 
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parameterize the background for mt,c with a function that behaves hke phase space near the 
endpoint (the ARGUS shape JlSl). In all fits, the signal and background yields are the only 
free parameters. 




-0.25 -0.125 0.125 0.25 
AE 



0.25-0.125 0.125 0.25 



FIG. 1. The AE distributions for B candidates (left) and B candidates (right). The solid line 
represents the sum of signal Gaussian (asymmetric for i^^vr^) and a linear background function 
determined from the fit. The signal component is represented by the dashed line. In the K^tt^ 
fits, the signal component is centered at —44 MeV due to vr mass assignment to the K track. A 
component centered at zero (dotted line) is also included to account for mis-identified B^ vr^vr" 
backgrounds. The K^n^ distributions are projections of a two-dimensional rribc and AE fit onto 
the AE distribution in the mbc signal region. 

Figure |l| shows the AE distributions for the B and B candidates. For K^n"^, the 
particle energy is calculated assuming a pion mass which shifts the signal peak to —44 MeV. 
A component centered at zero is also included to account for mis-identified B^ tt+tt" 
backgrounds. For the K^tt^ and Kgir^ distributions, the correct mass assignments are 
used. The results are listed in Table |. For the combined K'^tt asymmetry, the number 
of events is N{B) = N{K~n^) + ■yi^N{K~7i~^). In all cases, the measured asymmetry 
is consistent with zero. The results are also compared to the predictions in ||^ assuming 
03 = 90°. Other authors expect asymmetries that are less than 10%. 

The systematic error in the fitting procedure is determined by varying the fixed param- 
eters of the signal and background functions within their errors. The parameters are varied 
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TABLE I. Results of searches for DCPV in i? — > Kir decays. Listed are the number of signal 
events for each final state, the ^cp values with error, their 90% confidence intervals, and theoretical 
estimates from Ref. [5]. The first error is statistical; the second is systematic. The 90% statistical 
confidence intervals have been increased by the systematic error. In the K^tt^ final states, the 
asymmetry is corrected for the dilution {l — 2w), due to double mis-identification. For the combined 
K^TT asymmetry, the number of events are combined as N{B) = N{K~Tr^) + j^^N{K~7r~^). 



Mode 



N{B) 



N{B) 



i C.L. 



Theory [|] 



97 7+6.8 
^'•'-6.1 

16.511? 



5.6 



A"T(7r=^ + 7r°) 44.5 



+3.4 
-2.7 
+8.7 
-7.7 



h9.1 

-8.1 



44.2^ 



0.044" 



-0.186 +0.018 
-0.167 -0.021 



n nf;Q+0-222 +0.055 
-U.Ujy_Q j^gg _o,oi7 

-0.430 +0.020 
-0.343 -0.063 
-0.142 +0.017 
-0.126 -0.014 



0.098: 
0.003 



-0.25 
-0.40 
-0.53 
-0.22 



0.37 
0.36 
0.82 
0.25 



-0.19 
-0.18 
-0.01 



simultaneously for the conjugate final states and the deviation in ^cp is taken as the sys- 
tematic error. These deviations are on the order of ±0.01. We estimate the error due to 
other rare B decay backgrounds such as B ^ pn or K*'n' from the change in ^cp when the 
negative AE region is excluded from the fit. These backgrounds are shifted in AE by the 
mass of the missing particle, typically —140 MeV. This is the dominant systematic error in 
each mode: —0.017 for K^tt^, -1-0.05 for K^n^, and —0.06 for K^tt^. All deviations are 
then added in quadrature. 



> 3500 
2 3000 
S 2500 
w 2000 
I 1500 
1000 
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1.78 1.82 1.86 1.9 1.94 1.78 1.82 1.86 1.9 1.94 

FIG. 2. The K^tt^ invariant mass distributions for D meson candidates. The D peaks are 
fitted with a a = 5.5 MeV Gaussian and a second a = 13.6 MeV Gaussian containing 26% of the 
area. The solid curve is the sum the signal Gaussians and a linear background function (dashed 
line). The asymmetry in the yield is (1 it 6) x 10~^. 

We study the effects of detector-based asymmetries using a sample of high momentum 
D — > Ktc decays. Figure p| shows mass distributions for K~Ti^ and —>■ K^tt~ 

candidates. The daughter tracks are required to have |piab| > 1-5 GeV/c and the same 
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PID requirements as signal tracks. The D mass peak is fitted with a double Gaussian and 
a linear background function. We float all parameters in the fit to test for differences in 
resolution between the conjugate states. The signal and background shapes all agree within 
the errors for the conjugate states. We measure 19091 ± 172 candidates and 19055 ± 164 
candidates giving an asymmetry of (1±6) x 10"'^. We study this further using an inclusive 
sample of 2.2 x 10^ charged tracks (/i^) |]rB| required to have the same momentum and angular 
distributions as signal tracks. We measure an asymmetry of {N{h'^) — N{h~))/{N{h~^) + 
N{h~)) = (8.1 ±6.7) X 10~^, which is consistent with Monte Carlo expectations. We test for 
asymmetries in the PID performance by applying kaon and pion selection cuts to the charged 
hadron sample. The sample contains 1.4 x 10^ identified pions and 0.7 x 10^ identified kaons. 
We measure asymmetries of (A^(7r+) - A^(7r-))/(A^(7r+) + A^(7r")) = (1.6 ± 0.9) x 10"^ and 
{N{K+) - N{K-))/{N{K+) + N{K~)) = (-0.9 ± 1.1) x lO^^. These asymmetries result 
from a small charge dependence in the dE/dx measurement. Taking these numbers into 
account, we conservatively assign a 1% systematic error for detector-based asymmetries in 
all modes. This is added in quadrature to the fitting systematics discussed above, as well 
as the error in the dilution factor where appropriate. We add the systematic errors to the 
statistical 90% confidence intervals to obtain the final 90% confidence intervals listed in 
Table |. 

In conclusion, we have performed a search for DCPV in self tagging B Kit decays. All 
results are consistent with a null asymmetry and we conclude at the 90% confidence level that 
-0.25 < AcpiK^n"^) < 0.37, -0.40 < AcpiK^n^) < 0.36, and -0.53 < AcpiKln^) < 
0.82. For the combined K^n"^ and K^tt^ data samples we determine —0.22 < Acp{K^{7r'^ + 
TT^)) < 0.25. These results are consistent with the limits set in [§. The excellent performance 
of the Belle high momentum PID allows us to reduce significantly the dilution due to double 
mis-identification. This will become a crucial factor as the Belle data samples increase 
significantly. 
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Ministry of Education of Korea and the CHEP SRC program of the Korea Science and 
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No.2P03B 17017; the Ministry of Science and Technology of Russian Federation; the National 
Science Council and the Ministry of Education of Taiwan; the Japan- Taiwan Cooperative 
Program of the Interchange Association; and the U.S. Department of Energy. 
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